A&A manuscript no. 

(will be inserted by hand later) 



Your thesaurus codes are: 
11.03.5; 11.09.3; 11.05.2; 11.01.1 



ASTRONOMY 
AND 
ASTROPHYSICS 
1.2.2008 



OO 
ON 



> 

m 
O 
oo 
On 

Or 
i 

O 

c3 



13 



On the trend of [Mg/Fe] among giant elliptical galaxies 

Francesca Matteucci 1 ' 2 , RafFaella Ponzone 1,2 , Brad K. Gibson 3 

1 Dipartimento di Astronomia, Universita di Trieste, Via G.B. Tiepolo, 11, 34131 Trieste, Italy 

2 SISSA/ISAS, Via Beirut 2-4, 34014 Trieste, Italy 

3 Mt. Stromlo and Siding Spring Observatories, Australian National University, Weston Creek, Australia ACT 2611 
February 1, 2008 



Abstract. We revisit the problem of the flat slope of the 
Mg2 versus < Fe > relationship found for nuclei of ellipti- 
cal galaxies (Faber et al. 1992; Worthey et al. 1992; Carollo 
et al. 1993; Davies et al. 1993), indicating that the Mg/Fe 
ratio should increase with galactic luminosity and mass. 
We transform the abundance of Fe, as predicted by clas- 
sic wind models and alternative models for the chemical 
evolution of elliptical galaxies, into the metallicity indices 
Mg% and < Fe >, by means of the more recent index 
calibrations and show that none of the current models for 
the chemical evolution of elliptical galaxies is able to re- 
produce exactly the observed slope of the < Fe > versus 
Mg2 relation, although the existing spread in the data 
makes this comparison quite difficult. In other words, we 
can not clearly discriminate between models predicting a 
decrease (classic wind model) or an increase of such a ra- 
tio with galactic mass. The reason for this resides in the 
fact that the available observations show a large spread 
due mostly to the errors in the derivation of the < Fe > 
index. In our opinion this fact prevents us from drawing 
any firm conclusion on the behaviour of Mg and Fe in 
these galaxies. Moreover, as already shown by other au- 
thors, one should be careful in deriving trends in the real 
abundances just from the metallicity indices, since these 
latter depend also on other physical parameters than the 
metallicity. This is an important point since abundance 
ratios have been proven to represent strong constraints 
for galaxy formation mechanisms. 
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1. Introduction 

Elliptical galaxies do not show the presence of HII re- 
gions and it is not possible to resolve single stars in them 
in order to measure photospheric abundances. Therefore, 
most of the information on these objects is obtained from 
their integrated properties: abundances are derived either 
through colors or integrated spectra and in both cases the 
derived information is a complicated measure of metal- 
licity and age (the well known age-metallicity degener- 
acy). The most common metallicity indicators are Mgi 
and < Fe >, as originally defined in Faber et al. (1977; 
1985). Population synthesis techniques are adopted to an- 
alyze the integrated properties of ellipticals and to derive 
an estimate of their real abundances. Unfortunately, they 
contain several uncertainties residing either in incomplete 
knowledge of stellar evolution or in deficiencies in stel- 
lar libraries, as discussed in Chariot et al. (1996). In re- 
cent years more and more population synthesis models 
(Bruzual and Chariot, 1993; Buzzoni et al., 1992; Bres- 
san et al. , 1994; Gibson, 1997; Bressan et al., 1996; Gib- 
son and Matteucci, 1997; Tantalo et al., 1998) have ap- 
peared but the basic uncertainties still remain. In this pa- 
per we want to focus our attention about the comparison 
between theoretical model results and metallicity indica- 
tors. In this framework we will analyze the relationship 
between < Fe > and Mgi and its implications for the 
mechanism of galaxy formation. 

Several authors (Faber et al. 1992; Worthey et al. 
1992; Carollo et al. 1993; Davies et al. 1993; Carollo 
and Danziger 1994), from comparison of the observed in- 
dices with synthetic indices, concluded that the average 
[< Mg/Fe >]* in giant ellipticals must be larger than the 
solar value. This result was also confirmed by the analysis 
of Weiss et al. (1995) who made use, for the first time, of 
stellar evolutionary tracks calculated under the assump- 
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tion of non-solar ratios. The same authors found that the 
< Fe > versus Mg 2 relation among nuclei of giant ellip- 
ticals is rather flat and flatter than within galaxies. From 
the flat behavior of < Fe > vs. Mg 2 the same authors 
inferred that the abundance of Mg should increase faster 
than the abundance of Fe among nuclei of giant ellipti- 
cals. This conclusion is at variance with the predictions of 
supernova-driven wind models of ellipticals (Arimoto and 
Yoshii, f 987; Matteucci and Tornambe f 987). In fact, Mat- 
teucci and Tornambe (1987) showed that, in the frame- 
work of the classic wind model for ellipticals, the [Mg/Fe] 
ratio is a decreasing function of the galactic mass and lu- 
minosity. The reason for this behavior is clear: if Fe is 
mostly produced by the supernovae of type la, as it seems 
to be the case in our Galaxy (Greggio and Renzini 1983a; 
Matteucci and Greggio 1986), whereas Mg is mostly orig- 
inating from supernovae of type II, then the iron produc- 
tion is delayed relative to that of Mg and its abundance 
should be larger in more massive galaxies which develop 
a wind later than the less massive ones. All of this is valid 
under the assumption that after the onset of a galactic 
wind star formation should stop or should be negligible, 
which is a reasonable assumption for elliptical galaxies. 
Faber et al. (1992) proposed alternative scenarios to the 
classic supernova driven wind model, as originally pro- 
posed by Larson (1974). They suggested three different 
scenarios all based on the assumption that Mg is produced 
by type II supernovae and Fe is mostly produced by type 
la supernovae: i) a selective loss of metals, ii) a variable 
initial mass function (IMF) and iii) different timescales for 
star formation. These hypotheses were discussed by Mat- 
teucci (1994), who tested them in the context of chem- 
ical evolution models. In the hypothesis of the different 
timescales for star formation Matteucci (1994) suggested 
that the more massive ellipticals might experience a much 
stronger and faster star formation than less massive ellipti- 
cals leading to a situation where the most massive objects 
are able to develop galactic winds before the less massive 
ones. She called this case "inverse wind model" . On the 
other hand, in the classic wind model of Larson (1974) 
the efficiency of star formation was the same for all galax- 
ies thus leading to the fact that the galactic wind in more 
massive systems occurs later than in less massive ones, due 
to their deeper potential well. In the models of Arimoto 
and Yoshii (1987) and Matteucci and Tornambe (1987) the 
efficiency of star formation was a decreasing function of 
galactic mass, based on the assumption that the timescale 
for star formation is proportional to the cloud-cloud col- 
lision timescale which, in turn, is proportional to the gas 
density. Therefore, since in this monolithic collapse pic- 
ture the gas density decreases with the galactic mass, the 



galactic wind was even more delayed for the most massive 
systems. Matteucci (1994) proposed, as an alternative, a 
star formation efficiency increasing with the galactic mass 
and she justified this assumption by imagining that giant 
elliptical galaxies, instead of forming through a monolithic 
collapse of a gas cloud, form by merging of gaseous proto- 
clouds. The merging process can, in fact, produce higher 
densities for increasing galactic mass and/or higher cloud- 
cloud collision velocities resulting in a faster star forma- 
tion process. In such a model the galactic wind occurs ear- 
lier in massive than in smaller ellipticals thus producing 
the expected trend of an increasing [Mg/Fe] as a func- 
tion of galactic mass. Matteucci (1994) also showed that 
a variable IMF with the slope decreasing with increas- 
ing galactic mass and luminosity can produce the same 
effect without an inverse wind situation. The reason for 
that resides in the fact that a flatter IMF slope favors 
massive stars relative to low and intermediate masses, 
thus favoring Mg production over Fe production. How- 
ever, Matteucci (1994) could not translate the predicted 
abundances into Mg 2 and < Fe > since there were no 
available calibrations for [Fe/H] versus < Fe > but only 
calibrations for [Fe/H] vs. Mg 2 . Therefore she did not 
compare the predicted abundances with observations. 

Recently, calibrations for the iron index have become 
available (Borges et al., 1995; Tantalo et al., 1998) and 
therefore in this paper we revisit the whole problem of 
inferring trends on the real abundances by metallicity in- 
dices and we discuss the influence of the calibration re- 
lationships, which allow us to pass from indices to abun- 
dances, and we show that the inferred trend of Mg/Fe with 
galactic mass is not so clear when interpreted in terms 
of real abundances, thus warning us from drawing any 
firm conclusion on galaxy formation processes just on the 
basis of the observed behavior of < Fe > versus Mg 2 . 
indices. The reason for that resides partly in the large 
spread present in the observational data and partly in the 
fact that metallicity indices depend not only on the abun- 
dances of single elements but also on the ages and on the 
metallicity distribution (Tantalo et al. 1998) of the differ- 
ent stellar populations present in elliptical galaxies. 

In Section 2 we will discuss the chemical evolution 
model, in Section 3 we will define the average abundances 
of a composite stellar population, in Section 4 we will de- 
scribe the model results and transform the predicted abun- 
dances into indices by means of the most recent metallicity 
calibrations. Finally in Section 5 some conclusions will be 
drawn. 

2. The chemical evolution model 
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2.1. Basic equations 

The adopted model of galactic evolution is that outlined 
by Matteucci and Gibson (1995) (hereafter MG95), where 
extensive descriptions and references can be found. The 
evolution of the abundances of several chemical species 
(H, He, C, N, O, Ne, Mg, Si and Fc) in the gas is followed, 
taking into account detailed nucleosynthesis from stars of 
all masses and SNe of types la, lb, and II. We assume 
that ellipticals can be considered initially as homogeneous 
spheres of gas with luminous mass in the range 10 9 — > 
1 x 10 12 M . A single zone interstellar medium (ISM) 
with instantaneous mixing of gas is assumed throughout. 
The adopted age for all galaxy models is to, — 15 Gyr. 
The fundamental equations can be written 



dGjjt) 
dt 



-i>{t)Xi{t) - 
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1p(t - Tm)Qmi(t - r m )<f)(m)dm 
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J M Btn 

+ / ip{t - T m )Q mi (t - T m )cj>(m)dm, (1) 
J Mbm 

where Gi(t) = p gas (t)Xi(t) / p(Q) is the volume gas den- 
sity in the form of an element i normalized to the initial 
total volume gas density. The quantity Xi(t) represents 
the abundance by mass of an element i and by definition 
the summation over all the elements present in the gas 
mixture is equal to unity. 

The various integrals in equation (1) represent the 
rates at which SNe (I and II) as well as single low and 
intermediate mass stars and single massive stars restore 
their processed and unprocessed material into the ISM 
(for a detailed description of the integrals see MG95; Mat- 
teucci and Greggio 1986). We only remind here that the 
quantity Q m i represents the fraction of a star of mass m 
which is restored in to the ISM in the form of an element 
i and therefore contains the nucleosynthesis prescriptions 
that we assume to be the same as in MG95. 

The star formation rate tp(t) is given by: 



m = ^gas(i)MO). 



(2) 



i.e. normalized to the initial total volume density. ip(t) 
is assumed to drop to zero at the onset of the galac- 
tic wind. The quantity v is expressed in units of Gyr -1 
and represents the efficiency of star formation, namely 
the inverse of the time scale of star formation. The values 
adopted here for v are the same as in Matteucci (1992) and 
Arimoto and Yoshii (1987). In one case we have adopted 



the prescription for the inverse wind model of Mat- 
teucci (1994). The difference between the two cases is that 
in Matteucci (1994) the classic wind model assumes 
that the efficiency of star formation decreases with increas- 
ing total galactic mass, as in Arimoto and Yoshii (1987), 
whereas in the inverse wind model the efficiency of star 
formation increases as the total galactic mass increases 
(similar in spirit to the SFR efficiency parametrization of 
Tinsley and Larson 1979) leading to a situation in which 
more massive galaxies experience a galactic wind before 
the less massive ones. 

The quantity r m represents the lifetime of a star of 
mass m, and is taken from Padovani and Matteucci (1993). 



2.2. Galactic winds 

For gas to be expelled from a galaxy the following con- 
dition should be satisfied: the thermal energy of the gas 
heated by SN explosions should exceed the binding energy 
of the gas (Larson 1974). At this point the gas present in 
the galaxy is swept away and the subsequent evolution 
is determined only by the amount of matter and energy 
which is restored to the ISM by the dying stellar gener- 
ations. In particular, only low mass stars contribute to 
this evolutionary phase and, among the SNe, only SNe of 
type la (i.e. those SNe events whose progenitors have the 
longest lifetimes). 

Therefore, in order to evaluate the time for the on- 
set of a galactic wind we need to know the energy input 
from SNe and the binding energy of the gas as a function 
of time. The total thermal energy of the gas at time t, 
-E-thsN > an d the binding energy of the gas in presence of a 
diffuse halo of dark matter are calculated as described in 
Matteucci (1992). In particular, E t h SN (t) is calculated by 
assuming that ~ 70% of the initial blast wave energy is 
transferred into the ISM as thermal energy by a SN rem- 
nant, if the time elapsed from the SN explosion is shorter 
than a SN remnant cooling time (Cox 1972). The percent- 
age of transferred energy then decreases as a power law in 
time oc i~ - 62 for times larger than the cooling time. 

This is the same formulation used by Arimoto and 
Yoshii (1987), Matteucci and Tornambe (1987) and MG95 
with the exception that we consider a more realistic cool- 
ing time (expressed in years) which takes into account that 
the gas density is changing with time and that part of the 
interstellar gas is in the form of He: 



t c = 5.7 • 10 4 [2.99 • lCr 3 (M gas (i)/10 12 M Q ) . 

(M iuro /1O 12 M )- L6 t ' 53 years. 



(3) 
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Other formulations of the energy input from SNe can 
be found in Gibson (1997) and Gibson and Matteucci 
(1997). 

The energetic input from stellar winds in massive stars 
is ignored, since for normal ellipticals it is negligible com- 
pared to the SN thermal energy contribution, as shown by 
Gibson (1994). 

The galactic wind may last only for few 10 8 years or 
continue until the present time depending crucially on the 
assumptions about the thermal energy of the gas and the 
potential energy of the gas. Unfortunately, none of these 
quantities is well known. The time for the occurrence of 
the galactic wind, £gw, either increases with the galactic 
mass as a consequence of the potential well incresing with 
galactic mass together with the efficiency of star forma- 
tion decreasing with galactic mass (classic wind model), 
or decreases with the galactic mass if the efficiency of star 
formation is strongly increasing with galactic mass (in- 
verse wind model), as can be seen in Table 3, as will be 
discussed in Section 4. 

3. The average metallicity of a composite stellar 
population 

In order to compare model results with observations we 
should first calculate the average abundances of Mg and Fe 
of the composite stellar population of the given galaxies. 
The average metallicity or abundance in general which 
should be compared with the indices should be averaged 
on the visual light, namely: 



< Xi > L = 



E»j n t3 XiL Vj 
11, . ":,'-y 



(4) 



where n,j is the number of stars in the abundance 
interval Xi and luminosity Ly i . On the other hand, the 
real average abundance should be the mass-averaged one, 
namely: 



<X t > 



M- 



i r Sl 

S~Jo 



Xi{S)dS 



(5) 



where the subscript 1 refers to the specific time t\ (the 
present time) and Si is the total mass of stars ever born. 
Here we will use eq. (5) in order to compare model re- 
sults with indices and the reason is that we want to 
compare our results with previous ones (Matteucci 1994; 
MG95) and because for giant ellipticals the difference be- 
tween the mass-averaged metallicity and the luminosity- 
averaged one is negligible (Yoshii and Arimoto, 1987; Gib- 
son 1997). On the other hand, in smaller systems the dif- 
ference between the two abundances is not negligible due 



to the contribution to the light of low metallicity red gi- 
ants (Greggio 1996). We will show in Section 4 that in the 
cases studied here < Xi >l—< Xi >m- 

4. Model results 

In this section we will show the results of several models 
and we will convert, by means of the available calibrations, 
the average stellar abundances of Mg and Fe, [< Mg/H > 
]* and [< Fe/H >]*, as predicted for elliptical galaxies 
of different mass, into the metallicity indices Mg2 and 
< Fe >, respectively. 

Let us discuss first the metallicity calibrations. Re- 
lations linking the strength of metallicity indices to real 
abundances can be either empirical or theoretical. In the 
past few years several attempts have been made to cal- 
ibrate the strength of Mg 2 against [Fe/H] which has al- 
ways been considered as the measure of the "metallicity" 
in stars. It should be said that this is not entirely correct 
since we know that Mg does not evolve in loekstep with 
iron in the solar neighborhood nor in elliptical galaxies, 
due to the different timescales of production of these two 
elements. It would be much better to calibrate Mg 2 ver- 
sus [Mg/H] in order to avoid confusion. Calibrations of 
Mg2 versus < Fe > are from Mould (1978), Burstcin 
(1979), Peletier (1989), Buzzoni ct al. (1992), Worthey et 
al. (1992). In all of these calibrations the ratio [Mg/Fe] is 
assumed to be solar, at variance with the indication aris- 
ing from population synthesis models showing an over- 
abundance of Mg relative to iron in the nuclei of giant 
ellipticals (Faber et al. 1992; Worthey et al. 1992; Davies 
et al. 1993; Weiss et al. 1995). 

More recently, Barbuy (1994), Borges et al. (1995) and 
Tantalo et al. (1998) took into account non solar ratios of 
[Mg/Fe] in their calibrations. In addition, some of them 
(Borges et al. 1995; Tantalo et al. 1998) produced syn- 
thetic < Fe > indices thus allowing us to calibrate [Fc/H] 
also against < Fe >. This allows us to transform [Fe/H] 
in to Mg 2 and < Fe >, although many uncertainties arc 
involved in this exercise, mainly because, in this way, the 
derivations of Mg 2 and < Fe > are not independent. 

We run several models, in particular: Model I, which 
is the classic wind model, as described in MG95, with 
a Salpetcr (1955) IMF (namely an IMF with power in- 
dex x=1.35 over a mass range 0.1 < M/M Q < 100); 
Model II, which is the classic wind model with the Ari- 
moto and Yoshii (1987) IMF (namely an IMF with power 
index x=0.95 over the same mass range of the Salpeter 
one); Model III, which is the equivalent of the inverse 
wind model, as described in Matteucci (1994) with the 
Arimoto and Yoshii (1987) IMF; Model IV which is the 
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Table 1. Model I-classic wind, x=1.35 



Mi um (M e ) 




R eff (kpc) 


taw (Gyr) 




(Me) 




10 9 


19.0 


0.5 


0.131 


0.658 


■10 9 




10 10 


14.6 


1.0 


0.286 


0.791 


■ 10 10 




10 11 


11.2 


3.0 


0.514 


0.974 


■10 11 




10 12 


8.6 


10.0 


0.955 


0.885 


•10 12 




Table 2. Model II- classic wind, x=0.95 








M lum (M e ) 


v(Gyr- x ) 


R eff (kpc) 


taw (Gyr) 


M fin 


(M s ) 




10 9 


19.0 


0.5 


0.069 


0.246 


■10 9 




10 10 


14.6 


1.0 


0.117 


0.388 


■ 10 10 




10 11 


11.2 


3.0 


0.403 


0.668 


■10 11 




10 12 


8.6 


10.0 


0.660 


0.913 


■10 12 




Table 3. Model III- inverse wind, x= 


0.95 








Mi um {M Q ) 


v(Gyr- x ) 


R e ff (kpc) 


taw (Gyr) 


Mfi n 


(Me) 




10 9 


2. 


0.5 


1.200 


0.502 


■10 9 




10 10 


5. 


1.0 


0.897 


0.605 


■ 10 10 




10 11 


11. 


3.0 


0.408 


0.665 


■10 11 




10 12 


20. 


10.0 


0.205 


0.841 


•10 12 




Table 4. Model IV-classic wind, variable IMF 








M lum {M Q ) 




Reff (kpc) 


taw (Gyr) 


Mf in 


(Mq) 


IMF slope 


10 9 


19.0 


0.5 


0.163 


0.936 


■ 10 9 


2.0 


10 10 


14.6 


1.0 


0.237 


0.945 


■ 10 10 


1.7 


10 11 


11.2 


3.0 


0.514 


0.974 


■10 11 


1.35 


10 12 


8.6 


10.0 


0.660 


0.912 


•10 12 


0.95 


Table 5. Model V- classic wind, time variable IMF 








Mi um (M e ) 


v(Gyr- x ) 


Reff (kpc) 


taw (Gyr) 


Mj ln 


(M e ) 




10 9 


19.0 


0.5 


0.073 


0.276 


■10 9 




10 10 


14.6 


1.0 


0.223 


0.764 


■ 10 10 




10 11 


11.2 


3.0 


0.601 


0.960 


■10 11 




10 12 


8.6 


10.0 


2.205 


0.995 


•10 12 
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Table 6. Model VI- inverse wind, x=0.8 



.}{...} 



Mi um (M e ) 


v(Gyr- v ) 


Reff (kpc) 


taw (Gyr) 


M fm (M Q ) 


10 9 


4.0 


0.5 


0.22 


0.320 ■ 10 9 


10 10 


18.0 


1.0 


0.110 


0.430 • 10 10 


10 11 


32.0 


3.0 


0.090 


0.550 ■ 10 11 


10 12 


75.0 


10.0 


0.030 


0.450 • 10 12 



equivalent of the model with variable IMF, as described in 
Matteucci (1994), which assumes that ellipticals of smaller 
mass have a steeper IMF than the more massive ones. It 
is worth noting that the slope of the IMF is kept constant 
inside a galaxy. In particular, we vary the slope of the IMF 
from the Salpctcr one to the Arimoto and Yoshii one pass- 
ing from a galaxy with initial luminous mass of 1O 11 M 
to a galaxy with 10 12 M Q . This particular assumption can 
reproduce the observed tilt of the fundamental plane seen 
edge-on, namely the increase of M/L versus L as observed 
by Bender et al. (1992). 

Model V assumes a time-variable IMF as suggested 
by Padoan et al. (1997) and will be discussed in a forth- 
coming paper. In this formulation the IMF slope varies 
as a function of gas density and gas velocity dispersion, 
favoring the formation of massive stars at early epochs. 

Model VI assumes a constant IMF with a slope x=0.8 
and a star formation efficiency which varies more strongly 
with the luminous mass than in Model III. The slope 
x=0.8 is the limiting slope that we can accept to ob- 
tain a realistic M/Lb ratio for ellipticals, as discussed in 
Padovani and Matteucci (1993). 

The model parameters are described in Tables 1-6 
where we list the luminous masses in column 1, the star 
formation efficiency (in units of Gyr^ 1 ) in column 2, the 
effective radius (in units of Kpc) in column 3, the time 
for the occurrence of the galactic wind (in Gyr) in column 
4 and the final galactic luminous mass in column 5. For 
model IV is shown also the slope of the IMF in column 6. 

We then calculate the average [< Fe/H >]* and 
[< Mg/Fe >]* for the stellar component of ellipticals by 
using eq. (6) and finally we transform these abundances 
to observed Mg 2 and < Fe > line indices. In Tables 7- 
12 we show the results for different models and for the 
calibration of Tantalo et al. (1998). In particular, in Ta- 
bles 7-12 we show the luminous mass in column 1, the 
[< Fe/H >]* in the second column, the [< Mg/Fe >]* in 
column 3 and in column 4 and 5 the Mg 2 and the < Fe > 
indices, respectively. As already said, only the calibrations 
of Tantalo et al. (1998) and Borges et al. (1995) allow us 



to transform [Fc/H] into < Fe > and therefore allow us to 
compare model results with the data showing the behav- 
ior of Mg 2 vs. < Fe > among nuclei of giant ellipticals. In 
particular, starting from the synthetic indices of Tantalo 
et al. (1998) calculated for a fixed [Mg/Fe] and a fixed 
[Fe/H] we derived calibration relationships of the type: 

Mg 2 = f([Fe/H],[Mg/Fe]) (6) 

<Fe>=g([Fe/H},[Mg/Fe}) (7) 

which allowed us to derive the indices for any [Fe/H] and 
[Mg/Fe]. The calibrations we have adopted are: 

Mg 2 = 0.233 + 0.217[Mg/Fe] + 
+ (0.153 + 0.120[Mg/Fe]) ■ [Fe/H] (8) 



< Fe >= 3.078 + 0.341[M 5 /Fe] + 

(1.654 - 0.307[Mg/Fe\) ■ [Fe/H] (9) 

In Fig. 1 we show the mctallicity indices obtained by 
means of the already mentioned calibrations compared 
with the data (Gonzalez 1993; Worthey et al. 1992; Car- 
olio and Danziger 1994a, b). As one can easily see the data 
present a large spread, mostly due to the uncertainties in 
deriving the < Fe > indices. In particular, in Fig la we 
show the observed and predicted behavior of < Fe > vs. 
Mg 2 when Model I is adopted. The bestfit to the data 
implies the following relation, < Fe >= 3.94Mp 2 + 1-83, 
and is indicated in the figure. However, the spread in the 
data is large and this prevents us from drawing strong 
conclusions about a possible trend. The dotted lines in 
Fig. la represent the predictions of Model I obtained by 
means of the calibrations described before and adopting 
the same [Mg/Fe] ratio as predicted by the models, as one 
can see in Table 7. The agreement with the observed trend 
is not so good, showing that the slope of the predicted re- 
lation is steeper than that shown by the data and that 
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Table 7. Model I- classic wind, x=1.35 



log Mi um (M Q ) [< Fe/H >], (dex) [< Mg/Fe >]« 


(dex) 


Mp 2 (mag) 


< Fe > (A) 


9. -0.395 0.262 

10. -0.190 0.224 

11. -0.095 0.187 

12. -0.007 0.130 




0.217 
0.247 
0.257 
0.260 


2.545 
2.853 
2.990 
3.112 


Table 8. Model II- classic wind, x=0.95 


logM lum (M Q ) [< Fe/H >]» (dex) [< Mg/Fe >]« 


(dex) 


Mg 2 (mag) 


< Fe >(A) 


9. -0.833 0.397 

10. -0.034 0.329 

11. 0.388 0.288 

12. 0.502 0.264 




0.152 
0.298 
0.368 
0.383 


1.937 
3.137 
3.783 
3.957 


Table 9. Model III- inverse wind, x=0.95 


logM lum (M Q ) [< Fe/H >]» (dex) [< Mg/Fe >]» 


(dex) 


Mg 2 (mag) 


< Fe >(A) 


9. 0.197 0.262 

10. 0.334 0.274 

11. 0.384 0.290 

12. 0.317 0.308 




0.326 
0.355 
0.368 
0.360 


3.477 
3.700 
3.777 
3.678 


Table 10. Model IV- classic wind, variable IMF 


logMi um (M Q ) [< Fe/H >]» (dex) [< Mg/Fe >]» 


(dex) 


Mg 2 (mag) 


<Fe>(A) 


9. -1.450 0.149 

10. -0.834 0.169 

11. -0.095 0.187 

12. 0.502 0.264 




0.002 
0.125 
0.257 
0.368 


0.800 
1.800 
2.990 
3.957 


Table 11. Model V- classic wind, time variable IMF 


logMi um (M Q ) [< Fe/H >]» (dex) [< Mg/Fe >]» 


(dex) 


Mg 2 (mag) 


< Fe >(A) 


9. -0.550 0.339 

10. 0.115 0.253 

11. 0.267 0.136 

12. 0.291 -0.058 




0.200 
0.309 
0.308 
0.263 


2.341 
3.344 
3.555 
3.545 
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Table 12. Model VI- inverse wind, x=0.8 



.}{...} 



logM lum (M e ) 


[< Fe/H >], (dex) 


[< Mg/Fe >]» (dex) 


Mg 2 (mag) 


<Fe>(A) 


9. 


-0.099 


0.340 


0.287 


3.040 


10. 


-0.091 


0.365 


0.294 


3.064 


11. 


-0.023 


0.368 


0.308 


3.168 


12. 


-0.005 


0.390 


0.316 


3.188 



Fig. 1. a) Predicted and observed metallicity indices. The dot- 
ted line and open squares represent the < Fe > versus Mg 2 
predicted by Model I for galaxies of different masses, as indi- 
cated in Table 7. The dashed line and stars represent the real 
abundances of Fe and Mg as predicted by Model I for galaxies 
of different masses and arbitrarily translated in the < Fe > 
vs. Mg 2 diagram. The error bars referring to the different data 
samples are also shown, b) Predicted and observed Mg 2 versus 
mass relation. The predictions are from Model I. 



Fig. 2. a) The same as Fig. la but relative to the predictions of 
Model II. b) The same as Fig. lb but relative to the predictions 
of Model II. 

the predicted Mg 2 values do not cover the entire range in 
Mg 2 . This is mostly due to the assumed IMF since Model 
II, which assumes a flatter IMF, predicts values for Mg 2 
which cover the whole range (see Fig. 2a). 

In Fig. lb we show the predicted and observed mass- 
metallicity (Mg 2 ) relationship. The data are from Car- 
olio et al. (1993). The best-fit to these data indicate 
Mg 2 = 0.02ZogM tot + 0.08, where M tot is the total galactic 
mass (dark+luminous). The classic wind model recovers 
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the slope of the Mg 2 - mass relation, but with a zero- 
point offset of AMg 2 — 0.05 with respect to the observed 
distribution. 

On the other hand, the classic wind model with the 
Arimoto and Yoshii (1987) IMF (Model II), as shown 
in Fig. 2b, predicts a slope much steeper than the ob- 
served one, although it agrees better than Model I with 
the < Fe > vs. Mg 2 relation shown in Fig. 2a. It is worth 
noting that the Arimoto and Yoshii (1987) IMF well repro- 
duces the abundances in the intergalactic medium (MG95; 
Gibson 1997; Gibson and Matteucci 1997). It is worth not- 
ing that in figures la and 2a and in all the others we 
show also the relation between real abundances predicted 
by our models. The relation between [< Fe/H >]* and 
[< Mg /H >] * , arbitrarily translated in the plot of < Fe > 
versus Mg 2 , is indicated by the dashed lines. This is done 
only with the purpose of comparing the slope of the re- 
lation between real abundances with that of the relation 
between indices and they are very similar, indicating that 
the adopted calibration does not modify the predicted re- 
lation between Mg and Fe abundances. One of the main 
reasons for that is the adopted calibration which accounts 
for the right [< Mg/Fe >]* ratio for each galaxy. 

In Fig. 3a, b we show the predictions of the inverse 
wind model of Matteucci (1994) which predicts a stellar 
[< Mg/Fe >]* increasing with galactic mass. The slope 
of the < Fe > versus M g 2 relation is in better agreement 
than in the previous models, and the slope of the Mg 2 
vs. mass relation is also acceptable although the absolute 
values of the indices are too high. 

In Fig. 4a, b we show the results of Model IV with a 
variable IMF from galaxy to galaxy, which also predicts 
increasing [< Mg/Fe >]* ratios with galactic mass. The 
agreement with the < Fe > vs. Mg 2 data is marginally 
acceptable, but the slope of the mass-metallicity relation 
is too steep and the predicted absolute values of Mg 2 are 
too low. The low absolute values of Mg 2 are due to the 
fact that we used slopes steeper than the Salpcter (1955) 
one for the less massive galaxies and such slopes are not 
suitable for elliptical galaxies (sec MG95) since they pre- 
dict too low metallicities. However, other numerical ex- 
periments, where we used a variable IMF but with flatter 
slopes for each galactic mass (from x=1.4 in low mass 
galaxies to x=0.8 in high mass galaxies), have shown that 
there is a negligible difference in the predicted < Fe > 
vs. Mg 2 relation while the mass-metallicity relation gets 
worse. 

In Fig. 5a, b we show the predictions of Model V calcu- 
lated with the time-variable IMF as suggested by Padoan 
et al. (1997) and adapted to elliptical galaxies. The slope 
of this IMF is decreasing with time thus favoring massive 



Fig. 3. a) The same as Fig. la but for the predictions of Model 
III. b) The same as Fig. lb but for the predictions of Model 
III. 

stars at early epochs. A similar although more complex 
formulation of the Padoan ct al. (1997) IMF has been 
recently adopted by Chiosi ct al. (1998). The model be- 
haves like a classic wind model, in the sense that the 
galactic wind occurs first in small galaxies and later in 
the more massive ones. Concerning the predicted indices, 
Fig. 5a shows that the Mg 2 decreases for massive objects, 
due to the fact that the IMF in these galaxies is less bi- 
ased towards massive stars than in smaller systems. This 
is, in turn, due to the fact that the slope of the IMF is 
inversely proportional to the gas density which is lower in 
more massive objects. This model predicts a sort of bi- 
modal behavior for the Mg 2 vs. mass relation and it does 
not fit the data better than the other models. 

Finally, in Fig. 6a, b the predictions of Model VI are 
shown. At variance with all the previous models, Model 
VI seems to reproduce well the observed slope of the 
< Fe > vs. Mg 2 relation as well as the Mg 2 - mass re- 
lation. The main problem with this model is the fact that 
the predicted ranges of Mg 2 and < Fe > are too nar- 
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row compared to the observations, especially the range of 
< Fe >. Another potential problem of this model is the 
predicted Mj L b ratio which is ~ 30 for each galaxy mass. 
This is a high value for ellipticals unless one believes in a 
Hubble constant H a = 100Kms~ 1 Mpc~ 1 , as discussed in 
Padovani and Matteucci (1993). 



lead even to different slopes for the Mg 2 -logM relation- 
ship. 



Fig. 5. a) The same as Fig. la but for the predictions of Model 
V. b) The same as Fig. lb but for the predictions of Model V. 



Fig. 4. a) The same as Fig. la but for the predictions of Model 
IV. b) The same as Fig. lb but for the predictions of Model 
IV. 

In Figs. 7 and 8 we show the plot of the mass- 
metallicity (Mg 2 ) relation as predicted by Model I and 
Model II, obtained under different assumptions about the 
calibrating formula. As one can see, some of the calibra- 
tions give similar results such as those of Worthey ct al. 
(1992), Casuso et al. (1996) and Buzzoni ct al. (1992). 
These calibrations have in common the use of a solar ratio 
for [Mg/Fe] ([Mg/Fe]=0). On the other hand, the values of 
the indices obtained by using the calibrations of Barbuy 
(1994) and Tantalo et al. (1998) which adopt non-solar 
ratios, differ from the others and between themselves. It 
is worth noting that the use of different calibrations may 



One criticism that could in principle be moved to the 
results discussed before is that we adopted mass-averaged 
metallicities and not luminosity-averaged metallicities, as 
it should be the case. In Fig. 9 we show the indices ob- 
tained by the luminosity- and mass-averaged metallicities 
calculated for the results of Model II, when the calibra- 
tion of Tantalo et al. (1998) is applied. The luminosity- 
averaged metallicities, computed with the photometric 
model of Gibson (1997), are systematically slightly lower 
than the others and the difference is larger for smaller 
galaxies, as expected. However, the slope is the same in the 
two cases, showing that the use of mass-averaged metal- 
licity for this kind of analysis is quite justified. 

5. Discussion and Conclusions 

In this paper we have discussed the relation between 
mctallicity indices, such as Mg 2 and < Fe >, and total 
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Fig. 8. Mg2-niass relationships as predicted by Model II by us- 
ing different calibrations, as indicated in the figure. The best-fit 
to the data is also shown. 



Fig. 6. a) The same as Fig. la but for the predictions of Model 
VI. b) The same as Fig. lb but for the predictions of Model 
VI. 



Fig. 7. Mt/2-mass relationships as predicted by Model I by 
using different calibrations, as indicated in Fig. 8. The best-fit 
to the data is also shown. 



mass in nuclei of ellipticals and its implications in terms 
of models of formation and evolution of elliptical galaxies. 

In order to do that we have transformed the average 
abundance of Fe in the composite stellar population of 
the galaxy, as predicted by different models of chemical 
evolution, into Mg 2 and < Fe > indices by means of the 
available calibrations. 

We have shown the results of classic wind models for 
ellipticals, such as those discussed by Arimoto and Yoshii 
(1987) and Matteucci and Tornambe (1987), as well as the 
results of models with variable IMF from galaxy to galaxy 
and with galactic winds occurring first in the more massive 
systems, implying that these systems are older than the 
less massive ones. We have found that it is not possible to 
establish clearly which kind of model should be preferred, 
first of all because of the large spread present in the data. 

Moreover, little difference is found in the predicted in- 
dices of models which predict a [< Mg/Fe >]» either 
increasing or decreasing with galactic mass, although the 
data seem to suggest an increase of this ratio with galactic 
mass larger than predicted by any of the models. 

On this basis, the classic wind model can not be con- 
sidered worse than the other models. Actually, the classic 
wind model with a flat constant IMF seems to be the 
only one which can reproduce the whole range of the ob- 
served indices. However, if we isolate the data from Gon- 
zalez (1993) and do not consider the others, then in order 
to reproduce the flat slope of the < Fe > versus Mg 2 
relation, as given by the best-fit of the data, one should 
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Fig. 9. Predicted and observed < Fe > versus Mg^. The pre- 
dicted indices, relative to Model II, are calculated by averaging 
on the mass (stars) and on the visual luminosities (squares) . 



assume that Fe among the nuclei of ellipticals is almost 
constant whereas Mg increases from less massive to more 
massive nuclei. This is not achieved by any of the mod- 
els presented here since it would require quite "ad hoc" 
assumptions especially concerning the type la SNe. 

From the numerical experiments performed in this pa- 
per we can say that a model which explains at the same 
time the mass-metallicity and the iron-magnesium rela- 
tion requires an inverse wind situation, with a strong in- 
crease of the star formation efficiency with galactic mass 
(i.e. Model VI), rather than a variable IMF from galaxy 
to galaxy, and an IMF with a slope x = 0.8. However, a 
model of this type is not able to reproduce the observed 
ranges of < Fe > and Mg 2 - We have also calculated mod- 
els where amount and concentration of dark matter in- 
creases, compatibly with the formulation of the potential 
energy of the gas, with decreasing galactic luminous mass 
(see Persic et al. 1996), with the net result of obtaining 
an "inverse wind" situation. The results are very similar 
to those of Model III. Therefore, to obtain a better agree- 
ment with observations one should invoke also in this case 
an increase of the star formation efficiency with galactic 
mass. This would certainly flatten the < Fe > vs. Mgi 
relation but it would further shrink the ranges of the pre- 
dicted indices. In fact, both an increasing star formation 
efficiency and a decreasing importance of dark matter with 
increasing luminous galactic mass can be viable solutions 



to achieve the situation of more massive ellipticals being 
older than less massive ones. 

In conclusion, it is quite important to establish the 
value of [Mg/Fe] from the observational point of view 
since abundance ratios, such as [Mg/Fe], represent an im- 
portant diagnostic to infer ages in galaxies, due to the 
different timescales for the Mg and Fe production. Gen- 
erally, a high [Mg/Fe] is interpreted as a young age and 
the upper limit for the age is given by the time at which 
the chemical enrichment from type la SNe starts to be- 
come important. This timescale depends not only on the 
assumed progenitors of type la SNe but also on the star 
formation history of the galaxy considered (see Matteucci 
1997) and for giant elliptical galaxies this timescale is of 
the order of tsNeia ~ 3 — 4 • 10 8 years and in any case it 
can not be larger than 1 Gyr also for smaller systems. This 
is at variance with what stated by Kodama and Arimoto 
(1997) who claim, on the basis of results concerning our 
Galaxy (Yoshii et al. 1996), that this timescale is of the 
order of 1.5-2.5 Gyr. This is indeed true for our Galaxy 
where the star formation history has been quite different 
than in ellipticals and it had been already pointed out in 
Grcggio and Renzini (1983b) and in Matteucci and Greg- 
gio (1986). This is a quite important point, both for the 
galactic chemical enrichment and for the predictions about 
SN rates at high rcdshift. 

Therefore, an enhanced [Mg/Fe] indicates that the pro- 
cess of galaxy formation must have been very fast thus fa- 
voring a monolithic collapse scenario rather than a merg- 
ing scenario. In this framework, a [Mg/Fe] ratio higher in 
more massive ellipticals than in less massive ones could be 
interpreted as due to their faster formation and evolution 
(see Matteucci 1994; Bressan et al. 1996). 

An independent way of estimating the ages of ellipti- 
cals, where for ages we intend the time elapsed from the 
last star formation event, is to study the Hp index. This 
index is, in fact, related to the age of the dominant stellar 
population, since it gives a measure of the turn-off color 
and metallicity. It can therefore be used to solve the age- 
metallicity degeneracy. Bressan et al. (1996), by analyz- 
ing the H/3 and other physical parameters in the sample 
of ellipticals observed by Gonzalez (1993), concluded that 
massive galaxies should have stopped forming stars before 
less massive ones, in agreement with the results of the in- 
verse wind model discussed here. Finally, we would like 
to point out that both models with a Salpeter IMF and a 
variable IMF have a potential problem in reproducing high 
[q/Fc] ratios in the intracluster medium (ICM), as shown 
by their low average [< a/Fe >]* ratios (see Tables 7-12). 
Therefore, in agreement with MG95 and Gibson and Mat- 
teucci (1997) we conclude that a flat IMF is required to 
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explain the high [a/Fe] ratios, as found by ASCA obser- 
vations (Mushotzsky 1994). 
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